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ABSTRACT 

We present a study of the star formation histories of the Lupus I, III, and IV clouds using the Herschel 70-500 fim maps obtained by the Herschel 
Gould Belt Survey Key Project. By combining the new Herschel data with the existing Spitzer catalog we obtained an unprecedented census of 
prestellar sources and young stellar objects in the Lupus clouds, which allowed us to study the overall star formation rate (SFR) and efficiency 
(SFE). The high SFE of Lupus III, its decreasing SFR, and its large number of pre-main sequence stars with respect to proto- and prestellar sources, 
suggest that Lupus III is the most evolved cloud, and after having experienced a major star formation event in the past, is now approaching the end 
of its current star-forming cycle. Lupus I is currently undergoing a large star formation event, apparent by the increasing SFR, the large number 
of prestellar objects with respect to more evolved objects, and the high percentage of material at high extinction (e.g., above A v « 8 mag). Also 
Lupus IV has an increasing SFR; however, the relative number of prestellar sources is much lower, suggesting that its star formation has not yet 
reached its peak. 
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1. Introduction 

In the current paradigm of low-mass star formation (SF), a grav- 
itationally bound prestellar core will evolve into a young stel- 
lar object (YSO), passing through several phases, usually de- 
fined by Classes r epresenting increa s ing stages of evolution : 0, 
I, II, and III (see I Andre etal1l200Ql: lLada & Wilkingll984l for 
the definitions), before becoming a main- sequence star. While 

■ the later stages of low-mass SF are largely understood, less is 
known about the earlier stages (including the prestellar cores and 

, the Class objects) due to a lack of sensitivity and resolution 
at far-infrared t o submm waveleng ths. The Herschel Gould Belt 
Survey (HGBS, Andr e et al.l201 0). carried out with the Herschel 
Space Observatory (Pilbr attet al.ll2010l) . aims at studying these 
early stages of SF in nearby molecular clouds forming the so- 
called Gould Belt (IComer6netal.lll992h . 

Located at a distance betwe en 150pc (Lupl and IV) and 
200 pc (LupIII: IComeronl 120081) . the Lupus clouds (I, III, IV) 
are among the nearest star-forming regions in the Gould Belt. 
The large angular extent of the Lupus clouds across the sky 
(334°</<352°, 5°<Z?<25°) corresponds to a physical extent of 
50pcx55pc at a distance of 150 pc, similar to the distance 
range a mong the Lu pus I, III, and IV, clouds (50 pc). Previous 
Spitzer (jMerm et al. 2008, hereafter M08) and molecular-line 
dBenedettini et al.l l2012h studies of Lupus I, III, and IV have 
found that the three clouds seem to be at different stages of evo- 
lution: Lupus I is thought to be the youngest cloud, Lupus IV is 
a little more evolved, and Lupus III is the most evolved cloud. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



The Lupus I, III, and IV clouds were mapped, as a part of 
the HGBS, at five wavelengths from 70/mi to 500 /mi, covering 
the range where the spectral energy distribution (SED) of cold 
dust emission from prestellar sources and envelopes of Class 0/1 
objects (protostars), is likely to peak. Therefore, the Herschel 
data are crucial for detecting these objects and determining their 
physical parameters. The combination of the M08 catalog, con- 
taining mostly Class II/III pre-main sequence (PMS) stars, with 
the prestellar cores and Class 0/1 sources detected by Herschel 
allows for a much more complete view of the ongoing SF in 
the Lupus clouds than previously possible. Here, we present the 
first-look letter of the SF history of the Lupus clouds. A detailed 
analysis of the Herschel data on Lupus and the resulting cata- 
log of identified objects will be presented in a forthcoming first- 
generation paper (Benedettini et al. in prep.) and will be publicly 
available at the HBGS websiteQ. 

2. Herschel observations and data reduction 

The Lupus maps were obtained between January 2010 
and January 2011 by photometric observations with the 
Photodetector Array Camera and Spectrometer (PACS; 
Pog litsch et al.l |2010|) and Spectral and Photometric Imaging 
Receiver (SPIRE; [Griffin et all 1201 Oh in parallel mode using a 
scanning speed of 60"s _1 . Map sizes are 2°x2?3 for Lupus I 
and l?5xl?l for Lupus III, covering a similar region as the 
Spitzer observations for both clouds. Lupus IV was imaged in 
two maps, 2°xl?3 and l?3xl?3, covering the Spitzer map, as 
well as a new region, never mapped before at wavelength s of 
160-500/mi. The data were exported from HIPE v7.0 dOtu 
1201 Oh at level 0.5, and processed with the ROMAGAL data 

1 http://gouldbelt-herschel.cea.fr/archives 
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reduction pipeline (Traficante et al. 2011; Piazz o et al . 2012). 
The maps were astrometrically aligned with the 7 um MIPS 
image s from the Spitzer 'cores2disks' (c2d) survey (lEvans et al.l 
12003b , which in turn have been aligned with 2MASS data, based 
on a number of point sources observed in both 70/im maps 
and yielding a n astrometric preci sion of ~2"5. Absolute flux 
calibration (see Ber nard et al.ll2010l) was found to be better than 
20% by comparing the Herschel data with the Planck and IRAS 
data in the same area. The resulting maps have beam sizes of 9", 
12", 18", 25", 36", for 70/rni, 160 /mi, 250/im, 350/mi, and 
500 yum, respectively. The Lupus I observations were affected 
by stray Moonlight, visible as a bright band in the declination 
direction of the map. Fortunately, this did not affect our compact 
source fluxes, since the background is removed in the procedure. 

For computing column density and temperature maps, the 
70-350//m maps were convolved to the 500 yum resolution and 
rebinned to the same pixel size (1 1"5). The pixel-by-pixel modi- 
fied black-body fits were then performed on the regridded maps, 
excluding the 70 jim map since this emission might not be trac- 
ing the cold dust exactly. For the modified black- body fitting 
we a ssumed a dust opacity of Ki.3mm=0.004 cm 2 g _1 dHildebrandl 
Il983l) . a grain emi ssivity parameter /3=2 and the mean molecu- 
lar weis ht yu=2.8 (Kauffmann et al. 2008). The resulting column 
density map (Fig. lA.lb was used to define the cloud area: the 
emission within the Ay=2mag (assumi ng the column den sity to 
A v relation of Mt 2 =9.4x10 20 A v cm" 2 . iBohlin et af]ll978h con- 
tour was considered cloud emission and the integrated column 
densities in this contour (Table [B agree within 20% with the 
masses from the c2d extinction mapfl ( Chap man et al.l 12009b 
within the same area. 

Compact source detection and extraction were performed 
with CuTEx ( Molinar i et al.ll20lH) in each of the five Herschel 
maps sep a rately, using a 3cr SNR detection limit. Following 
Elia et al.l (120101) . sources across the five bands were associ- 
ated according to their positions. We adopted a conservative ap- 
proach and removed sources that had a displaced counterpart in 
the 350yum or 500 yum bands by more than half of the FWHP 
from the source center common to the other bands, which can 
introduce a large uncertainty in the measured fluxes at the longer 
wavelengths. In particular, for the prestellar cores this could re- 
sult in finding false objects by an overestimation of the mass 
(see next section). Sources with detection in fewer than three 
bands longward of A=10 jim were discarded. We then fit a sin- 
gle temperature modified black-body function to the SEDs of 
the extracted sources, again excluding the 70 yum flux, to derive 
dust temperatures. Source masses were then derived from the 
optically thin part of the modified black-body spectrum, using 
the same dust properties as above. The reduced chi- squared was 
used to remove sources with badly fitted SEDs, which would 
yield uncertain masses and temperatures otherwise, and hence 
could have a faulty prestellar core classification. Objects iden- 
tified as non-Lupus mem bers based on their proper motions 
dLopez Marti et al.l 1201 ll) . in total 16 objects, were removed. 
We cross-checked our sample with the SIMBAD and the c2d 
database to remove known extragalactic sources (in total 13 
galaxies were found). Seven off-cloud unresolved objects with- 
out a 70yum counterpart were classified as candidate galaxies 
and removed. Finally, to understand whether some of the off- 
cloud candidate protostars might have been misclassified (see 
Fig. lA.ll and Sect. [4]), we estimated the galaxy coun t complete- 
ness. Using the number counts from lClements et ID 420 1 0), we 
find that for Lupus I and III we start lacking galaxies with flux 



2 http://data.spitzer.caltech.edu/popular/c2d/2007 1 101 _enhanced_vl/ 



levels of ~350mJy at 250/mi, while for Lupus IV the limit is 
700 mJy. 



3. Results 



Table 1. Properties of the Lupus clouds 



Cloud 


Lupus I 


Lupus III 


Lupus IV 


distance 61 (pc) 


150 


200 


150 


coverage (degree 2 )^ 


4.6 


1.6 


4.3 


cloud area>A v =2mag (pc 2 ) 


16.4 


7.7 


7.1 


cloud mass>Ay=2mag (M©) 


830 


570 


500 


cloud mass>Av=8 mag (M©) 


145 


65 


60 


Mot = Af Y SO + Afp r estellar 


52 


113 


37 


N YS0 0/I/II/IIF 


1/10/10/4 


1/10/50/42 


3/9/11/4 


prestellar 


27 


10 


10 


N unbound 


68 


8 


63 



Notes. {a) See lComeroiil d2008h {b) SPIRE and PACS overlapping area. 
(c) Number of objects per Class 0, 1, II, and III are based on the Herschel 
data and on the M08 catalog. 



Figure IA.1I shows that the Lupus clouds are overall elon- 
gated in shape and fragmented into smaller clumps, in which 
most of the pre- and protostellar objects are located. A prestel- 
lar core is defined as a gravitationally bound dense core (size 
< 0.05 pc) without an internal luminosit y source that will sup- 
posedly form stars in the future (see e.g., Wa rd-Thompson etaLl 
[l99l|d i Francesco et al. I l2007h . We used the ratio of the core 
mass to its critical Bonnor-Ebert mass to distinguish be- 
tween prestell ar cores and unbound objects that may or may not 
form stars (cf. IKonyves et al.ll20Tol) . wher e M is the cor e mass, 
M BE =2ARc 2 s /G the Bonnor-Ebert mass dBonnoril 19561) . c s the 
sound speed (c s = ^/Jcb T l\imu), R the radius, and T the core tem- 
perature. For the core mass, we used the mass obtained from 
the modified black-body fitting and computed the M B e using the 
fitted modified black-body temperature and the deconvolved ra- 
dius at 250//m. Cores with ^ > 1.0 are gravitationally bound 
dBonnorl 19561) and were classified as prestellar, while cores with 
a lower ratio were considered unbound objects. An example of 
prestellar core SED is given in Fig.[T] 

A protostar is defined as a core that is internally heated by an 
embedded YSO and is still accreting material from its envelope. 
The pre sence of a YSO is asc ertained by a 70 /mi and/or 24 jim 
source (IKonyves et al.l 12010). Among the protostellar objects , 
the submm to bolometric luminosity ratio ( Andre et al.ll2000l) . 

^s2b = ^tubm - m ^bo\ Afim deriv ed & om the Herschel, M08, and 
WISE data ( WrighI eTaTll201Ql) . was used to identify the proto- 
stars as Class or Class I objects. Protostars with L S 2t>>3% were 
classified as candidate Class 0, while protostars with a lower 
ratio were identified as candi date Class I (a more c onservative 
Class definition than used by lBontemps et al.ll2010l) . Examples 
of a Class and a Class I SED are shown in Fig.Q] The classifi- 
cation of more evolved objects, i.e., the PMS stars (or Class II/III 
objects), was taken from M08 and SIMBAD. M08 also identified 
flat SED objects tha t trace the transition between Class I and II 
dGreene et al.ll 19941) . Since we cannot distinguish between Class 
I and F sources with the Herschel data, we considered them as 
Class I objects. 

With the conservative source selection (Sect. 2) and the clas- 
sification described above, we found 47 candidate prestellar 
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Fig. 2. Five-band images of three prestellar cores (contour levels are peakflux*0.99, 0.95, 0.90, 0.80, 0.70, 0.60, and 0.50) and 
a Class object in Lupus I. The maps are centered on the prestellar core/Class object and are ordered by right ascension. The 
remaining prestellar and Class sources of all three Lupus clouds are given in Fig. lA.2[ 
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Fig. 1. Spectral energy distributions of a prestellar core, Class 0, 
and a Class I source in Lupus I and the modified black-body fit to 
the cold dust emission (blue line). The 3.4-24 yum and 1300 yum 
data points were taken from M08. 



cores (see Fig.[2l continued in Fig. IA.21 for all the five-band im- 
ages). For the entire Lupus complex, the prestellar cores have 
median masses of 0.25 M Q within a range of 0. 1-3.0 M , and 
median temperatures of 9 K within a range of 7-14 K. The mean 
prestellar source si ze is 0.02 ± 0.01 pc. Five Class sources were 
found (Figs. [2l and IA.2I) , of which most were previously classi- 
fied as Class I objects (M08). The on ly previously known Class 
is Lup MM3 dTachihara et alJl2007h . The Class temperatures 
are between 9.5 K and 15 K. Their envelope masses ranged from 
0.08 to 0.6 M Q . For two Class sources, we compared the dust 
temperatures obtained from the He rschel data to the kinetic tem- 
peratures obtained from ammonia ( Benede ttini et al. 2012), and 
found that they agree within 20%. We found 14 new Class I can- 



didates in the Lupus clouds, 1 1 of which were outside the field 
covered by Spitzer c2d. 

Table Q] summarizes the properties of the Lupus clouds and 
the objects found therein based on our results and the M08 cata- 
log. The overlap between the Herschel data, including the YSOs 
that were only detected at A=l 0/im and 160/zm, with the M08 
catalog was quite large, especially for the younger Class I ob- 
jects. We detected 61% of their Class I objects, 46% of their 
Class II objects, and 10% of their Class III objects in the area 
covered by both surveys. For the few undetected Class I objects, 
we place an upper limit of 66 mJy at 250 yum, and assuming a 
point source morphology, this implies a mass limit of 0.002 M© 
at a dust temperature of 15 K. 



4. Star formation history 



iPalla & Stahierl {2000) posit that the star formation rate (SFR) 
has been increasing in the Lupus clouds over the past 4Myr 
using number counts of PMS stars, whose ages they estimated 
through comparison with theoretical evolutionary PMS tracks. 
Our sample contains both prestellar cores and Class 0/1 objects, 
and by merging with the M08 catalog we can, therefore, estimate 
the recent SFR behavior better than in previous studies. 

One can deduce the relative number of objects N of a cer- 
tain class with respect to a reference class expected for a con- 
stant SFR as N=N Y QfXr/r Y Qf by assuming the classes' lifetimes 
r and r re f. Since the Class II objects are expected to be com- 
plete (M08 and SIMBAD), we chose this as the reference clas s 
with the lifetime estimate of 2+1 Myr (see lEvans et al . 2009). 
For the other classes, we assume 0.5 Myr for the prestellar 
lifetime (lEnoch et al.l[2008h . 0.05 Myr for the Class lifetime 
([Froebrich et al.1 120061) . and 0.84 Myr for the Class I lifetime 
(lEvans et al.ll2009l) . Using these lifetimes, we derived the ex- 
pected numbers of objects in each Class (relative to the num- 
ber of Class II objects). Figure [3] shows the ratio 77 of the 
observed-to-expected source numbers for the prestellar (Tj pres ), 
Class (?7cio), and I (77a 1) objects in each cloud. In this figure, 
we also plot the 7701 without the possible galaxy contaminated 
sources (see Fig. IA.lb based on the estimated completeness limit. 
Clearly, this possible contamination does not influence our SFR 
conclusions. 

For both Lupus I and IV we find more prestellar, Class 0, and 
Class I objects than predicted for a constant SFR. We are there- 
fore witnessing an increasing SFR over the past 0.5-1.5 Myr. 
The 77's, particularly the 77 pres , is much larger in Lupus I than in 
Lupus IV, suggesting that Lupus I is undergoing a star formation 
event, while the SFR of Lupus IV might increase even more in 
the future. In Lupus III there are fewer prestellar, Class 0, and 
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Fig. 3. Ratio of observed-to-expected (for a constant SFR) 
source numbers, 77, per Class. Solid lines represent the numbers 
from Table [TJ while dashed lines show Class I without the pos- 
sible galaxy contamination. 



Class I objects than expected for constant SFR. The star forma- 
tion in Lupus III has decelerated over the past 2 Myr. 

The SFR analysis depends heavily on the lifetimes used in 
the estimation. Performing the inverse analysis, assuming a con- 
stant SFR, and deriving the lifetimes (with respect to the number 
of Class II objects) would lead to lifetimes longer by a factor ~10 
for the prestellar sources, ~4 for the Class objects, and ~2 for 
the Class I objects in Lupus I. Finally, the results of this analy- 
sis do not depend on the source-finding algorit hm, CuTEx, since 
analy sis of the sources found by getsources (iMen'shchikov et all 
l2012h yielded similar SFR trends. 

We calculated the star formation efficiency (SFE), which is 
the ratio of the total mass in YSOs (Class - III) M s , assuming 
0.2 M Q as the average Lupus PMS star mass (M08), to the total 
mass of the cloud plus YSOs: M Ms +M • With this formulation, 
we find a similar SFE of ~1% in Lupus I and IV, but a much 
higher SFE of 3.5%in Lupus III. The latter may be expected from 
the high YSO-to-prestellar objects ratio (Table [T]). The different 
distance of Lupus III, 200 pc, does not influence these results, 
because putting Lupus III at 150 pc magnifies its SFE to 6.1%. 

One can study the star formation histories by interpreting 
the SFR behavior and SFE of the clouds as different evolution- 
ary states within a star formation cycle. We propose that the high 
SFE of Lupus III, its decreasing SFR, and large number of PMS 
stars with respect to proto- and prestellar sources suggest that 
Lupus III is the most evolved cloud, which after having expe- 
rienced a major star formation event, is now approaching the 
end of its current star-forming cycle. The propertie s of Lupus III 
seem to be similar to those of Chameleon I, where Bell oche et alJ 
(2011) claim to see the end of star formation based on prestel- 
lar cores found in 870 /mi LABOCA data. On the other hand, 
Lupus I is currently undergoing a large star formation event, ap- 
parent by the increasing SFR and the large number of prestel- 
lar objects with respect to more evolved sources. Also Lupus IV 
has an increasing SFR; however, the relative number of prestel- 
lar sources is much lower, suggesting that its star formation has 
not yet reached its peak and that Lupus IV is at an earlier stage 
of evolution than Lupus I. The contrast between the increasing 
SFR in Lupus I and the decreasing SFR in Lupus III is possi- 
bly reminiscent of the contrast between LI 688 and LI 689 in 
Ophiuchus. There, L 1688 is an active star-forming cloud with 



many prestellar cores, while L 1689 contains only a few prestel- 
lar cores, even though these clouds have similar CO properties. 
However, while iNutter et al] (120061) explain this difference by an 
external trigger from a nearby OB association, we claim that the 
diverse SFRs and SFEs in Lupus result from the clouds being in 
different states of their star formation cycle. 

In Table Q] we also list the fraction of cloud mass above 
the SF threshold of Ay^8mag found in recent studies (e.g., 
iHeiderman et al.ll2010t lAndre et al.ll2010h . This number, which 
is independent of source selection and identification, provides 
an estimate of the percentage of cloud mass directly available 
for SF, and it is expected to s cale a s the probability of finding 
prestellar cores. Spezzi et al.1 (1201 11) were the first to use this 
method on the extinction maps of the Lupus clouds. We ap- 
plied this method to the Herschel data and found that Lupus 
I has a much higher percentage of mass above the Ay thresh- 
old than the Lupus III and IV clouds, supporting the idea of a 
star formation event in this cloud. Lupus IV has a only slightly 
higher p ercentage of dense m aterial than Lupus III, as was also 
noted by Spe zzi et all (1201 ll) . supporting the similar number of 
prestellar cores found in these two clouds. The fraction of cloud 
mass above this Ay threshold correlates well with the number of 
prestellar cores found with Herschel, thereby strengthening our 
conclusions. 

Acknowledgements. The authors thank the anonymous referee for her/his com- 
ments. SPIRE has been developed by a consortium of institutes led by Cardiff 
Univ. (UK) and including: Univ. Lethbridge (Canada); NAOC (China); CEA, 
LAM (France); IFSI, Univ. Padua (Italy); I AC (Spain); Stockholm Observatory 
(Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, Univ. Sussex 
(UK); and Caltech, JPL, NHSC, Univ. Colorado (USA). This development has 
been supported by national funding agencies: CSA (Canada); NAOC (China); 
CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); 
STFC, UKSA (UK); and NASA (USA). PACS has been developed by a con- 
sortium of institutes led by MPE (Germany) and including UVIE (Austria); 
KU Leuven, CSL, IMEC (Belgium); CEA, LAM (France); MPIA (Germany); 
INAF-IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); IAC (Spain). This devel- 
opment has been supported by the funding agencies BMVIT (Austria), ESA- 
PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI/INAF (Italy), 
and CICYT/MCYT (Spain). KLJR, ES, DE, MP, and DP are funded by an ASI 
fellowship under contract numbers 1/005/1 1/0 and 1/038/08/0. NLJC is supported 
by the Belgian Federal Science Policy Office via ESAs PRODEX Program. 

References 

Andre, P., Men'shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L 102 
Andre, P., Ward-Thompson, D., & Barsony, M. 2000, PPIV, 59 
Belloche, A., Schuller, E, Parise, B., et al. 2011, A&A, 527, A145 
Benedettini, M., Pezzuto, S., Burton, M. G., et al. 2012, MNRAS, 419, 238 
Bernard, J.-R, Paradis, D., Marshall, D. J., et al. 2010, A&A, 518, L88 
Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 
Bonnor, W. B. 1956, MNRAS, 116, 351 

Bontemps, S., Andre, P., Konyves, V., et al. 2010, A&A, 518, L85 

Chapman, N. L., Mundy, L. G., Lai, S.-R, & Evans, N. J. 2009, ApJ, 690, 496 

Clements, D. L., Rigby, E., Maddox, S., et al. 2010, A&A, 518, L8 

Comeron, F. 2008, Handbook of Star Forming Regions II, ed. Reipurth, B., 295 

Comeron, E, Torra, J., & Gomez, A. E. 1992, Ap&SS, 187, 187 

di Francesco, J., Evans, II, N. J., Caselli, P., et al. 2007, PPV, 17 

Elia, D., Schisano, E., Molinari, S., et al. 2010, A&A, 518, L97 

Enoch, M. L., Evans, II, N. J., Sargent, A. I., et al. 2008, ApJ, 684, 1240 

Evans, II, N. J., Allen, L. E., Blake, G. A., et al. 2003, PASP, 115, 965 

Evans, II, N. J., Dunham, M. M., J0rgensen, J. K., et al. 2009, ApJS, 181, 321 

Froebrich, D., Schmeja, S., Smith, M. D., et al. 2006, MNRAS, 368, 435 

Greene, T. P., Wilking, B. A., Andre, P., et al. 1994, ApJ, 434, 614 

Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3 

Heiderman, A. and Evans, II, N. J. and Allen, L. E., et al. 2010, ApJ, 723, 1019 

Hildebrand, R. H. 1983, QJRAS, 24, 267 

KaufFmann, J., Bertoldi, E, Bourke, T. L., et al. 2008, A&A, 487, 993 
Konyves, V., Andre, P., Men'shchikov, A., et al. 2010, A&A, 518, L106 
Lada, C. J. & Wilking, B. A. 1984, ApJ, 287, 610 

Lopez Marti, B., Jimenez-Esteban, E, & Solano, E. 2011, A&A, 529, A108 
Men'shchikov, A., Andre, P., Didelon, P., et al. 2012, A&A, 542, A81 



4 



K. L. J. Rygl et al.: Recent star formation in the Lupus clouds as seen by Herschel 



Merin, B., J0rgensen, J., Spezzi, L., et al. 2008, ApJS, 177, 551, M08 
Molinari, S., Schisano, E., Faustini, R, et al. 2011, A&A, 530, A133 
Nutter, D., Ward-Thompson, D., & Andre, P. 2006, MNRAS, 368, 1833 
Ott, S. 2010, in ASP Conf. Series, Vol. 434, 139 
Palla, F. & Stahler, S. W. 2000, ApJ, 540, 255 
Piazzo, L., Ikhenaode, D., Natoli, P., et al. 2012, in IEEE, Vol. 21, 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, LI 
Poglitsch, A., Waelkens, C, Geis, N., et al. 2010, A&A, 518, L2 
Spezzi, L., Vernazza, P., Merin, B., et al. 2011, ApJ, 730, 65 
Tachihara, K., Rengel, M., Nakajima, Y., et al. 2007, ApJ, 659, 1382 
Traficante, A., Calzoletti, L., Veneziani, M., et al. 2011, MNRAS, 416, 2932 
Ward-Thompson, D., Scott, P. R, Hills, R. E., et al. 1994, MNRAS, 268, 276 
Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 



Appendix A: Online Material 



Istituto di Astrofisica e Planetologia Spaziali (INAF-IAPS), via del 
Fosso del Cavaliere 100, 00133 Roma, Italy 
e-mail: kazi . rygl@inaf . it 

2 Laboratoire AIM Paris-Saclay, CEA/IRFU CNRS/INSU Universite 
Paris Diderot, 91191 Gif-sur-Yvette, France 

3 CESR, Observatoire Midi-Pyrnes (CNRS-UPS), Universite de 
Toulouse, BP 44346, 31028 Toulouse, France 

4 Rutherford Appleton Library, Chilton, Didcot, OX11 0NL, UK 

5 Department of Physics and Astronomy, Open University, Milton 
Keynes, UK 

6 University of Athens, Department of Astrophysics, Astronomy and 
Mechanics, Faculty of Physics, Panepistimiopolis, 15784 Zografos, 
Athens, Greece 

7 CNRS/INSU, Laboratoire d'Astrophysique de Bordeaux UMR 
5904, BP 89, 33271 Floirac, France 

8 Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 
3001 Leuven, Belgium 

9 Department of Physics and Astronomy, University of Victoria, PO 
Box 355, STN CSC, Victoria BC Canada, V8W 3P6 

10 National Research Council Canada, Herzberg Institute of 
Astrophysics, 5071 West Saanich Road, Victoria BC Canada, V9E 
2E7 

11 Canadian Institute for Theoretical Astrophysics (CITA), University 
of Toronto, 60 St. George Street, Toronto ON Canada, M5S 3H8 

12 European Southern Observatory, Karl-Schwarzschild-Strasse 2, 
87548 Garching bei Munchen, Germany 

13 INAF-Osservatorio Astrofisico di Arcetri, Large E. Fermi 5, 50125 
Firenze, Italy 

14 Jeremiah Horrocks Institute, University of Central Lancashire, PR1 
2HE, UK 



5 



K. L. J. Rygl et al.: Recent star formation in the Lupus clouds as seen by Herschel, Online Material p 1 




Right Ascension (J2000) 




4 i i I i i i i i I I -] q 1 

16 h 10 m 0Q s 
Right Ascension (J2000) 



Fig. A.l. H2 column density maps of Lupus I, III, and IV in units of 10 20 cm" 2 with Ay contours overplotted. For Lupus III and 
IV the contours are Ay = 2, 3, 6, and 9 mag, while for Lupus I the contours are 4, 6, and 9 mag (to avoid the stray Moonlight). 
The different classes of objects, from the Herschel data and from the M08 catalog, are indicated. The off-cloud candidate Class I 
sources, marked by an asterisk, are those we considered as possible galaxy contaminations (shown by a dashed line in Fig. [3]). The 
one off-cloud Class I in Lupus I without an asterisk has strong 24 jim emission, so is less likely to be a galaxy. Dashed contours 
mark the Spitzer MIPS coverage. 
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Fig. A.l. - Continued 
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Fig. A.2. Visual catalog of the prestellar cores (contour levels are peakflux*0.99, 0.95, 0.90, 0.80, 0.70, 0.60, and 0.50) and Class 
objects in the three Lupus clouds (I, III, IV), continued from Fig El The maps are centered on the prestellar core/Class object and 
are ordered per cloud by right ascension. 
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Fig. A.2. - Continued 



